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Abstract 



We report on calculations of the reduced sedimentation velocity U /Uq in 
homogenous suspensions of strongly and weakly charged colloidal spheres 
as a function of particle volume fraction (p. For dilute suspensions of 
strongly charged spheres at low salinity, U /Uq is well represented by the 
parametric form 1 — p<p a with a fractional exponent a = 1/3 and a param- 
eter p ~ 1.8, which is essentially independent from the macroion charge 
Z. This non-linear volume fraction dependence can be quantitatively un- 
derstood in terms of a model of effective hard spheres with ^-dependent 
diameter. For weakly charged spheres in a deionized solvent, we show that 
the exponent a can be equal to 1/2, if an expression for U/Uq given by 
Petsev and Denkov [J. Colloid Interface Sci. 149, 329 (1992)] is employed. 
We further show that the range of validity of this expression is limited to 
very small values of <j) and Z, which are probably not accessible in sedi- 
mentation experiments. The presented results might also hold for other 
systems like spherical proteins or ionic micelles. 



1 Introduction 

The sedimentation velocity U of interacting colloidal particles depends both on 
the indirect hydrodynamic interactions (HI) mediated by the suspending solvent, 
and on the microstructure of the suspension. In equilibrium, the latter is de- 
termined by direct potential forces arising for example from the steric repulsion 
between the particles and from the electrostatic repulsion of overlapping double 
layers. Different pair potentials u(r) lead to rather different microstructures, and 
this in turn strongly effects the sedimentation. It is well established by theo- 
ries and experiments that the sedimentation velocity of a dilute suspension of 
monodisperse hard spheres is given by 

¥r = 1-6.550 + O(0 2 ), (1) 

where <fi is the particle volume fraction, and Uq is the sedimentation velocity at 
infinite dilution [0-0. On the other hand, the long-ranged electrostatic repulsion 
occurring in suspensions of charged particles can give rise to a reduction in U, as 



compared to a hard sphere dispersion at the same volume fraction [§-11]. This 
decrease in U is mainly due to the cumulative backflow of displaced fluid, which 
becomes particularly effective because the probability of two or more charged 
particles coming close to each other is very small. Conversely, hard spheres are 
effectively attracted to each other at small interparticle distances r ~ 3a, where 
a is the radius of the spheres. This can be readily seen from the potential of 
mean force w(r) = —kBTlng(r), where g(r) is the radial distribution function 
of hard spheres, which shows its maximum value at contact distance r = 2a 
due to excluded volume effects. Closely spaced particles are mutually exposed 
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to the downflow of nearby fluid, dragged along with the sedimenting particles. 
Consequently, the retardation from backflow is reduced, whereas the influence of 
near-field HI is enhanced for hard sphere suspensions. 

While the reduction in U due to long-ranged electrostatic forces is known since 
many years, it was not realized until recently that, in particular for deionized 
charge-stabilized suspensions, U can be significantly smaller than Uq even for 
extremely small volume fractions as <fi ~ 10~ 4 . In fact, in the past the effects 
of HI have been frequently considered to be negligible for dilute suspensions of 
charged particles |12| , |l3l| . However, recent calculations have clearly demonstrated 
for these systems that HI is of importance at small <fi not only for sedimentation 
1^-[T6|, but also for short-time |17[) |18|1 and long-time |L9| collective diffusion, 



and for long-time self-diffusion [20|. Moreover, it was shown theoretically for 
dilute suspensions of charged colloids without added electrolyte that U follows a 
non-linear ^-dependence of the parametric form flH|[H)| , |TS[ 

^ = l-p^. (2) 

The numerically calculated coefficient p ~ 1.8 was found to be nearly indepen- 
dent of the (effective) particle charge Z, provided that Z is kept large enough to 
completely mask the hard core of the particles. Experimenal results of sedimenta- 
tion experiments on charged colloids agree favorable with the scaling-prediction 
of eq. (0) |I6|) pi]]. Similar non-linear volume fraction dependencies are found 
for the short-time translational and rotational diffusion coefficients of charged 
particles p2}-|2"4|. 



The same ^-dependence as found for salt-free fluid suspensions of charged 
colloids is known to be valid for the sedimentation velocity of dilute ordered 



arrays of fixed spheres poTPq], For such arrays, the coefficient p is determined 
analytically as p = 1.76 for a sc lattice, and as p = 1.79 for a fee or a bee 
lattice |25l,p7||. As discussed in detail in Refs. fl4l, |l5|,|l8| , the main peak position 
r m of the radial distribution function g(r) for highly charged particles in salt-free 
suspension scales (p 1 / 3 which is also typical of a crystalline solid. This, in 

fact, turns out to be relevant for eq. (fj) to be valid for strongly repelling particles 
like charged colloids and for ordered arrays of fixed spheres |14| , [I^ ,|18|)p^]. 

In the present article, we analyse the sedimentation velocity of charged col- 
loidal dispersions at low salinity both as a function of and of Z. Our numerical 
results for U/Uq with system parameters representing monodisperse modified 



PMMA particles investigated very recently in sedimentation experiments [fZI 
are well described by eq. @ if the effective macroion charge is chosen such that 
Z > 100. We physically explain the scaling relation eq. (0) in terms of an effec- 



tive hard sphere (EHS) model by using Wertheims analytical expression [29| for 
the static structure function of hard sphere dispersions obtained in Percus-Yevick 
approximation. In this specific form of the EHS model, also the numerically de- 
termined value p ~ 1.8 is recovered very accurately. We will further demonstrate 
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for weakly charged spheres that U/Uq — 1 displays a volume fraction dependence 
proportional to 1//2 , provided the suspension is completely deionized and is 
chosen small enough that the radial distribution function of the macroions can 
be approximated by its zero-density limit. This result, however, turns out to 
be valid only for extremely small values of <fi and Z, which we believe are not 
accessible in sedimentation experiments. 



2 Theory of sedimentation in charged colloids 

We start this section by summarizing the theoretical method used to calculate 
the reduced sedimentation velocity U/UqoI charged colloidal spheres. Our results 
for U/Uq are based on the effective macroion fluid model of charge-stabilized 
suspensions In this model, the effective pair potential u(r) between two 

charged colloidal particles consists of a hard-core part with radius a, and of a 
longer-ranged screened Coulomb potential u e i(r) for r > 2a, i.e. 

e ~n(r-2a) 

f3u el {r) = 2Ka . (3) 

Here, K is a dimensionless coupling parameter given by 

K = hL(^J)\ (4) 
2a \1 + Kaj ' K J 

where L B = (3e 2 /e is the so-called Bjerrum length, and (3 = (fc^T) -1 is the 
thermal energy. The suspending fluid is treated as a continuum without internal 
structure, only characterized by its dielectric constant e. Moreover, Z is the 
effective charge of a colloidal particle in units of the elementary charge e. The 
screening parameter k is given by 

k 2 = 4tcL b \n\Z\ + 2n s ] = k\ + k? 8 , (5) 

where n s is the number density of added 1-1-electrolyte, and n = 30/(47ra 3 ) is 
the number density of colloidal particles. Notice that k comprises a contribution 
k c due to counterions, which are assumed to be monovalent, and a second con- 
tribution k s arising from added electrolyte. Eq. (H) is a good approximation for 
u(r) even for strongly charged colloids and for values of na significantly larger 
than one, provided the effective charge number Z is regarded as an adjustable 
parameter [ |T5|j . 

On the time scales probed by dynamic light scattering and by sedimentation 
experiments, the effect of HI on the translational motion of the colloidal particles 
is embodied in the hydrodynamic diffusivity tensors DfAR N ); i,j = 1, . . . , N 



12, [13, [L5|, |3(| |3T|| . The short-hand notation R = (Ri, . . . , Rn) denotes the 

tt i 

ij ( 



configuration of the N spherical particles. Without HI, DfjR N ) = 5ijD l 



4 



where Dq = ksT '/ '(pirrja) is the Stokesian diffusion coefficient of a particle with 
radius a in a solvent with viscosity rj, and 1 denotes the unit tensor. 

The reduced short-time sedimentation velocity U/Uq of a macroscopically ho- 
mogeneous suspension of monodisperse colloidal spheres is then given by the 
zero-wavenumber limit 

— = lim#(g) (6) 



U () 



of the so-called hydrodynamic function H(q) |13, 15, 171- This function is defined 

as 

1 



H(q) 



ND> 



i, 



JV 

E 



R 



q e 



(7) 



with wavevector q of magnitude q and corresponding unit vector q = q/q. The 
brackets indicate an equilibrium ensemble average. H(q) can be regarded as 
a generalized (short-time) sedimentation coefficient of particles exposed to spa- 
tially periodic external forces aligned with q, and derived from a weak potential 
proportional to exp [— iq ■ r] pO], |3~^1 . 

In principle, one needs to distinguish the short-time sedimentation velocity 
U, defined through eq. (|6]), from the long-time sedimentation velocity, which 
is determined in conventional sedimentation experiments. In dilute suspensions, 
however, when the HI are well described as a sum of pairwise additive interac- 
tions, both quantities are identical, since then sedimentation does not perturb 



the equilibrium microstructure |SB|, . At larger volume fractions, the mi- 
crostructure becomes distorted from its equilibrium form, since n-body HI with 
n > 3 becomes important. This distortion leads to an additional change in 
the sedimentation velocity. Nevertheless, simulations of hard sphere suspensions 
show that the differences between the short-time and the long-time sedimentation 
velocities due to memory effects are rather small f35|,|3"6|. 

The assumption of pairwise additive HI is justified for the important case of 
dilute (typically (p < 0.1) charge-stabilized suspensions at sufficiently low ionic 
strength, as considered here. For such systems, the particles are kept apart from 
each other due to their strong electrostatic repulsion. Then the iV-body diffusivity 
tensors Dfj(R N ) are well approximated by the two-body tensors 



D i f\R N ) = 6 lJ 



N 



1 + 

i=i 



R, 



1 - 8ij)B(Ri - Rj) 



where the term I = i is excluded from the sum. This approximation has been 
verified in the case of translational and rotational self- diffusion by considering also 
the leading three-body contribution to HI ||22,2^|, and in case of the sedimentation 



velocity by comparison with a more elaborate method, known as the lowest order 
form of the 57-expansion [fl4, |37| , |38! . The two-body translational mobility tensors 
A(r) and B[r) are calculated by means of series expansions in powers of (a/r) 
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3^.|40[. We only quote the leading terms for further reference 

15 /a\ 4 



A(r) = -— {-) rr + 0( r - b ) (9) 
B(r) = U-)[l + rr] + U-X[l-3rr] + 0(r- 7 ), (10) 



4 \rj 2 \r 

with r = r/r. Using eqs. (0,H), H(q) is expressed in terms of an integral 

H(q) — 1 + n J dr g(r) [q ■ A{r) ■ q + q ■ B{r) ■ q cos(q ■ r)\ , (11) 

involving these mobility tensors together with the radial distribution function 
g(r). In this work, we use eqs. (f|ll) together with the series expansions of A(r) 
and B{r) for calculating U/U , by including contributions up to order (a/r) 20 . 

When terms only up to 0(r~ 4 ) in the series expansions of A(r) and B(r) are 
employed, U/Uq is given explicitly as 



U 



5 + 3 I" dx x(l - g(x)) + — I" dx 9 ^ 

J2 4 J2 



X 2 



(12) 



with x = r/a. 

For charge-stabilized suspensions, it is not necessary to account for many 
terms in the expansion of the two-body mobility tensors, since the integrals in 
eq. (|i~T| ) converge rapidly because g(r) is essentially zero at small interparticle 
distances. In contrast, many terms are needed for hard spheres. For example, 
using the zero-density form go(r) = B(r — 2a) for the radial distribution function 
of hard spheres, we obtain from eq. ( |12D that U/Uq = 1 — p<ft with p = 6.87. Here, 
Q(x) is the unit step function. By including terms only up to 0(r~ 3 ), the result 
is p = 5.0, as can be seen from eq. (|i~2|), when the last term on the right hand 
side, arising from the term proportional to (a/r) 4 in eq. ® is omitted. On the 
other hand, if terms up to 0(r~ 20 ) are considered, the result for p is improved to 
p = 6.54, which is close to the exact value p = 6.55, first obtained by Batchelor 
using tabulated numerical results for the near-field HI ||T|. 

To obtain U/Uq for dilute charge-stabilized suspensions, we calculate g(r) in 
the effective macroion fluid model by using, for simplicity, the well-established 
rescaled mean spherical approximation (RMS A) |15| . On the basis of the pairwise 
additivity approximation of the HI combined with the (a/r)-expansion of the mo- 
bility tensors, henceforth referred to as PA-scheme, U/Uq can then be calculated 
as a function of 6. 



3 Fractional density dependence of U/Uq 
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3.1 Strongly charged particles 

We show in the following that the exponent 1/3 in eq. (0) and the charge- 
independence of the parameter p, both found from numerical calculations, can 
be understood quantitatively in terms of a model of effective hard spheres (EHS), 
which can be treated analytically. 

The EHS model accounts for the most important feature of the radial distri- 
bution function g(r) of highly charged particles, namely the so-called correlation 
hole. For an illustration of that, consider fig. [II which shows a typical g{r) 
for a salt-free suspension of strongly charged particles at = 0.08. Due to the 
strong electrostatic repulsion between the particles, g(r) has a well developed first 
maximum, and a spherical region with nearly zero probability of finding another 
particle. This region is referred to as the correlation hole. For small 0, the cor- 



relation hole usually extends over several particle diameters [jl5| ,p4. Therefore, 
we can approximate the actual g{r) of the charge-stabilized system by the radial 
distribution function of an effective hard sphere (EHS) system with an effective 
radius aEHS > a and an effective volume fraction 4>ehs = 4>{aEHs/a) 3 - The 
EHS-radius aEHS accounts for the electrostatic repulsion between the particles 
and can be identified as aEHS = f m /2, where r m is the principal peak position of 
the actual g(r). 

Since the extension of the correlation hole is substantially larger than a for 
small volume fractions and large particle charges, it is then a good approximation 



to consider only the leading Oseen-term in B(r) (cf. eqs. (|9|,10jTTD) in calculating 



the sedimentation velocity. U/Uq is then well approximated by 



' , . *> , 

Uq a 2 Jo a 



1 + drrh(r) = l + ^H(s = 0). (13) 



Here, h(r) = g(r) — 1 is the total correlation function and 

_ roa 

H(s) = / drrh(r)e~ sr (14) 



denotes the Laplace transform of rh(r). Next, we approximate h(r) by the total 
correlation function h,EHs( r ] (Pehs) of the EHS model, evaluated at the effective 
volume fraction 4>ehs- With this approximation for h(r) used in eq. (|13|), we 
readily obtain the result U/Uq = 1 — p<p 1 ^ with a fractional exponent 1/3, and 
the parameter p given by 

P = 34> 2 ^ s H E hs(z = 0), (15) 

where Hehs( z ) is the Laplace transform of y hEHs{y',<f>EHs) with y = r/aEHS 
and z = saEHS- 

We notice that <Pehs and hence p are indeed independent of (f> and Z (> 100), 
provided that Qehs is identified with r m /2. This follows from the fact that, for 



7 



deionized suspensions of strongly charged particles where k c k s holds, r m 
coincides within 3% with the average geometrical distance r = a(30/(47r)) _1//3 of 
two particles. For an illustration of this fact consider the inset of fig. [I], which 
shows RMSA results for r m /f as a function of (p. Thus, the scaling relation 
r m « r oc </) -1 / 3 holds and this gives rise to the exponent 1/3 in eq. (fj): By 
identifying clehs with r, we obtain an effective volume fraction 4>ehs = tt/6 
independent of Z and </>. 

To obtain a numerical value of p according to eq. flT5]), we take now advantage 
of an analytic expression for Hehs( z ) given in the Percus-Yevick (PY) approx- 
imation By performing the zero-z limit, we obtain in PY approximation 
after a straightforward calculation the intermediate result 

tlEHSV> - UJ , r • (IbJ 

5 (1 + 2(pEHS> 

Substitution of eq. (|16"D into eq. (|T5| ) gives 

5 1 + 2 s 



Since 0^H5 = 7r/6, we obtain from eq. (|T7|) a value p = 1.76 remarkably close 
to p — 1.80 as determined from the parametric fit of our numerical PA-result for 
U/Uq (cf. following section). 

For an illustration of the replacement of the actual g(r) by the EHS-g(r) 
within the EHS model, see again fig. |], which shows besides a typical RMSA- 
g(r) of charge-stabilized particles the corresponding EHS-g(r) obtained in PY 
approximation for 4>ehs — 71 / 6. 

We finish this section with two remarks. First, we would like to stress again the 
close connection between the relation eq. (0), found for highly charged colloidal 
suspensions, and the corresponding result for ordered arrays of fixed spheres. For 
both systems, the 1//3 -scaling behavior of U/Uq is caused by a strong structural 
correlation of the particles, resulting in the scaling relation r m oc </> -1 / 3 (|14|,|l"5 



18] . Although Saffmann already mentioned the possibility for finding a scaling 



behavior as eq. (g) for highly correlated fluid suspensions ||26|| , the numerical 
method descibed here provides the first quantitative results on the sedimentation 
velocity of charged suspensions, leading to the predicted 1//3 -scaling in eq. (fj) 
plus the calculation of the prefactor p. It is further important to notice that the 
scaling relation r m oc 0" 1 / 3 is not valid if significantly large amounts of excess 
electrolyte are added to the suspension, since then the particle diameter becomes 
another physically relevant length scale besides the mean particle distance f. The 
0-dependence of U/Uq — 1 at small (f) changes then gradually with increasing n s 
from a 1 / 3 -dependence to the linear ^-dependence of eq. (]!]) characteristic for 
hard spheres [IB] . 
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3.2 Weakly charged particles 



It is also of interest to consider the opposite limiting case of dilute suspensions of 
weakly charged spheres in the framework of the effective macroion fluid model. 
This case was addressed first by Petsev and Denkov |il,|42|. Suppose Z is so 
small that 

(3u el (r = 2a) = K < 1, (18) 

i.e. the electrostatic repulsion can be treated as a small perturbation of the 
hard-core repulsion between the particles. For very small 0, g(r) can then be 
approximated as 

g{r) ~ <7o(r) [1 - (3u el {r)\ (19) 

with go{r) = G(r — 2a). In other words, g(r) is approximated by its zero-density 
limit g(r) ~ exp (— f3u(r)) = go{r) exp (— /3u e i(r)), linearized with respect to 
f3u e i{r). Substitution of eq. (|19|) in eqs. (|B|,|TTD gives two additive contribu- 
tions to U /Uq, the first one arising from the hard-core part go(r), and the second 
one from the electrostatic perturbation of g(r) in eq. ([19]). Explicitly |f41, 42 



U 



K 

6.55 + 6— 

Ka. 



(20) 



where the first term in the bracket is the result of eq. (|l|), arising from the hard- 
core part of g(r). As far as the electrostatic part is concerned, it was argued 
by Petsev and Denkov that the leading Oseen-term in B(r) should give the 
dominant contribution to U /Uq. Therefore only the Oseen-term is considered 
here, leading to the second term in the bracket of eq. (p0|). We will provide a 
critical discussion of this approximation in the following section. 

Next, we further simplify eq. (p0|) for the case of vanishing excess electrolyte 
(n s = 0), where na = A<p 1 / 2 with A = (3\Z\Lb/ 'a) 1 / 2 : . When and/or \Z\ are 
sufficiently small so that Acj) 1 ^ 2 <C 1, eq. (|20|) is further simplified to 



U 
Uo 



— ~ 1 



2a 



+ 0(<f>). 



(21) 



Thus, we expect the reduced sedimentation velocity of dilute deionized suspen- 
sions of weakly charged particles to scale like U/Uo = 1 —pcf) 1 ^ 2 , with a parameter 
p depending on the macroion charge Z and on the ratio L j e/(2a). 

Notice that the hard core contribution 6.550 to U/Uq in eq. (|20|) is omitted 

a/2 



against the term proportional 1/2 in proceeding from eq. ( p0|) to eq. (pl|) . 
Therefore, eqs. (^) and (|2l|) do not become equal to each other for Z —>■ 0. 
Whereas eq. fl20| ) reduces to the hard sphere result of eq. ([]]) for Z = 0, eq. (f2~T|) 
results in U = Uq for vanishing particle charge. Obviously, the range of validity of 
eq. pip is restricted to small values of Z, but still sufficiently large that the second 
term in eq. (pO| ) plays the dominant role compared to the hard core contribution 
6.550. We mention already here that due to the strong approximations made in 
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deriving eqs. (|20|J21| ), the range of validity of these two equations is restricted to 
extremely small values of </> and Z (cf. following section). 

For eqs. (|20jj2T|) to be valid, it is further assumed that the van der Waals 
attraction among the particles is negligibly small. Experimentally, this might 
be achieved by coating the particle surface with a thin layer of polymer chains 
and/or by matching the solvent refractive index to the refractive index of the 
colloidal particles. If attractions between the particles are present, not only eqs. 
(2"0",pT) cease to be valid, but also sedimentation velocities larger than for hard 
spheres may be measured at the same density |5|,|6|, [l^ , |4^ ,H- This enhanced sed- 
imentation rate is due to the enhanced probability, as compared to hard spheres, 
of two particles for being close together, thus leading to a reduced retardation 
from backflow p3| . 



4 Numerical results and discussion 

In this section, we present numerical results for the reduced sedimentation ve- 
locity U/Uq as a function of the volume fraction <fi and of the particle charge 
number Z. Our PA-scheme calculations of U/Uq account, if not stated differ- 
ently, for two-body contributions to the hydrodynamic mobility tensors up to 
order (a/r) 20 . The system parameters employed in our calculations are e = 2.183 
(corresponding to cis-decalin as dispersing fluid at temperature T = 293K), and 
particle radius a = 695nm, representing modified PMMA particles investigated 
very recently in sedimentation experiments We further use n s = 0, i.e. the 
ionic strength is essentially determined by the (monovalent) counterions, which 
counterbalance the charge of the colloidal particles. Obviously, the dielectric con- 
stant e used in the calculations presented here is rather small, leading to strong 
long-ranged repulsions between the particles even for small surface charges, which, 
most probably, are present in any polar organic solvent [^5 |. Furthermore, for 
these solvents, residual water can cause strong electric charging of the PMMA 
particles [^TJ. We therefore present calculations for both weakly and strongly 



charged particles. Our calculations for strongly charged particles recover the 
results obtained in Refs. JT^, [1^], [T8[] for different system parameters, demonstrat- 



ing clearly the independence of the scaling relation eq. (g) from specific system 
parameters, in particular from the solvent dielectric constant e. 



4.1 Strongly charged particles 

In fig. 0, we show the PA-scheme result for U/Uq, obtained by choosing an 
effective charge number Z = 150 large enough that the physical hard-core radius 
a of the particles constitutes no relevant physical length scale. This numerical 
result is perfectly fitted by the form l—p<j) a , with p = 1.80 and a = 0.34 ~ 1/3. As 
mentioned before, this is in remarkably good quantitative agreement with our EHS 
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model result obtained in PY approximation. As discussed in detail in Ref. [|16 
adding small amounts of excess electrolyte leads to a significant increase in U, 
and the 1//3 -scaling behaviour of U/Uq does not hold any more. For comparison, 
fig. H] includes also the reduced sedimentation velocity of uncharged hard spheres 
according to eq. (|I|). Evidently, the sedimentation velocity of charged particles 
decreases much faster with increasing <f> than one would expect for hard spheres 
at the same volume fraction. 

As shown in fig. [I], the RMSA-g(r) corresponding to the largest volume frac- 
tion = 0.08 considered in fig. |^ has well developed undulations, with the 
maximum value approximately located at the mean geometrical particle distance 
r. We only quote that the corresponding static structure factor S(q) has its prin- 
cipal peak height well below the range of values 2.8 — 3.1, where the system starts 



to freeze according to the empirical Hansen- Verlet rule [fig] . A few comments on 
the quality of the RMS A input for g(r), used in the PA-scheme calculations, 
is in order here. It is well known |L5[] that the RMSA underestimates to some 
extend the oscillations of g(r) and S(q) in case of strongly correlated systems 
of highly charged particles. In principle, we could use instead of the RMSA an 
alternative scheme like the Rogers- Young (RY) scheme ||47|1 . The RY sheme is 
quite accurate within the effective macrofluid model, but has the disadvantage 
of being numerically far more involved than the RMSA. Fortunately, the RMSA 
has been found to give nearly identical results as the RY-scheme, provided that 
a somewhat larger value of Z is used in the RMSA calculations ||15|| . We have 
argued before and will show in the following that the parameter p in eq. (|2]) is 
nearly independent of Z, typically as long as Z > 100. Consequently, using the 
RY-scheme instead of the RMSA for the same value of Z should lead to practi- 
cally identical results for U/Uq. We have verified this assertion by explicit RY 
calculations |4^,|49fl. 

In fig. |3|, we display PA-scheme results for U/Uq versus for various values 
of Z. For Z = 0, we recover the hard sphere result U /Uq = 1 — p<p with p = 6.54. 
An increase in Z leads to a strong reduction in U, with a gradual transition from 
the linear 0-dependence of U/Uq towards the non-linear form of eq. (|2]), with 
p ~ 1.80. This figure nicely illustrates that U/Uq becomes independent of Z for 
Z > 100. For these large values of Z, r m stays practically constant when Z is 



increased with fixed d> 48 



The effect on the PA-results for U/Uq caused by truncating the (a/r )-expansion 
of the two-particle mobilities after various terms of increasing order in (a/r), can 
be assessed from fig. |j. The solid line represents the full PA-scheme result where 
all two-body contributions up to 0(r~ 20 ) are accounted for. Nearly identical 
results for U/Uq are obtained, even at <fi = 0.08, when contributions only up to 
0(r~ A ) are considered. In contrast, the result for U/Uq obtained by accounting 
only for hydrodynamic contributions up to 0(r~ 3 ) shows clear deviations from 
the full PA-scheme result at volume fractions <fi > 0.01. As a conclusion, we 
can state that for dispersions of highly charged particles it is justified to use a 
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truncated far-field expansion of the mobility tensors with only the first terms 
being included, provided a good static input for g(r) is used. The finding that 
an expansion up to C(r -4 ) leads already to good results indicates further that 
hydrodynamic n-body contributions to U with n > 3 are indeed negligibly small 
in the considered 0-range. For practical purposes, it is therefore suitable to use 
the simple eq. ([Tj^) for calculating the reduced sedimentation velocity for highly 
charged suspensions. 

4.2 Weakly charged particles 

We discuss now the behaviour of U/Uq when Z is small. To be specific, consider 
first the value Z = 5. For this charge number, f3u e i(r) < (Lb/ (2a))Z 2 ~ 0.5. This 
implies that the microstructure is affected not only by the electrostatic repulsion 
but also by the physical hard core of the particles. In fig. |5], we have redrawn 
from fig. the PA-scheme result for U/Uq with Z = 5 and RMSA-input for g(r). 
This graph should be compared with the corresponding result for U/Uq obtained 
from the expression given in eq. (^) (henceforth referred to as PD- result). The 
PD-result for U/Uq is completely different from the corresponding PA-scheme 
result, in particular at larger where the PD-U/Uo even turns negative. The 
failure of the PD-expression in describing U/Uq arises from the fact that the 
approximations entering in its derivation, i.e. in particular the approximation 
eq. (|i~9l) for g(r), but also the truncation of the electrostatic contribution to 
U/Uq after the leading Oseen-term, are not justified even for a charge number 
as small as Z = 5. The consideration of more terms in the hydrodynamic pair 
mobility contributions up to 0(r~ 20 ), however, leads only to a small increase in 
U/Uq, as can be seen from fig. |5], which further includes the PA-scheme result 
for U/Uq with g(r) approximated by the linearized zero-density form of g(r), eq. 
(|19|). Another slight improvement in U/Uo is achieved when the non-linearized 
zero-density form g{r) = 6(r — 2a) exp (— (5u e i{r)) is used as input in the PA- 
scheme. In fact, the contact value f3u e i(2a) is not small enough for Z = 5 to fully 
justify a linearization in f3u e i(r). Therefore, the main reason for the failure of the 
PD-expression is due to the significant differences of the actual g(r) at Z = 5 
from its zero-density form already for ~ 0.01. 

To obtain further insight in the range of validity of the PD-result eq. (^), we 
display in fig. |6| results for U/Uq for an even smaller particle charge Z = 2 and 
volume fractions smaller than = 0.01. As seen from this figure and from fig. || 
the differences between the RMS A- PA-scheme result for U/Uq and the PA-scheme 
result for using g(r) = 6(r — 2a) exp(— (3u e i{r)) or g(r) = 0(r — 2a) [1 — j3u e i(r)] 
as static input become smaller with decreasing Z, but are still significant even for 
Z = 2 and ~ 10~ 3 . For an explanation of this finding, we refer to fig. |7|, where 
the radial distribution g(r) obtained from the RMS A and from the zero-density 
expression g(r) = @(r — 2a) exp(— f3u e i(r)) for Z = 2 and = 0.005 are shown. 
Obviously, there are still small differences between the the zero-density expression 
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for g(r) and the corresponding RMS A result even for these small values of Z and 
0. We show now that these small differences in g(r) cause the large differences 
in the results for U/Uq as illustrated in fig. |6]. For this purpose, we plot in the 
inset of fig. [5] the function r(l — g{r)) as obtained from the radial distribution 
functions shown in the main figure. As easily seen, the two curves for r(l —g(r)), 
corresponding to the two different inputs for g(r), are remarkably different from 
each other even up to large distances r. It is now crucial to notice that according 
to eq. fll^D it is essentially the function r(l — g(r)) and not g(r) itself which 
appears in the integrand when U/Uq is calculated including terms up to 0{r~ l ) 
in the series expansions in eqs. (|9|JT0|). Therefore, the small differences in the two 
g(r)'s considered in fig. [7|give rise to substancial differences in the corresponding 
sedimentation velocities shown in fig. |6|. We can therefore conclude that it is 
absolutely important to employ an accurate g(r)-input for calculating U/Uq in 
the PA-scheme. Due to the use of the (linearized) zero-density approximation of 
g(r), the expression given by Petsev and Denkov in eq. (p0|) does not predict the 
sedimentation velocity of weakly charged particles correctly. 

After having explored the high sensibility of U/Uq on the form of the g(r)- 
input used in the PA-scheme, it is now apparent that the second approximation 
made in deriving eq. (p0|), i.e. the omission of terms in the far-field expansions of 
the hydrodynamic mobilities of 0(r~ 3 ) in the electrostatic contribution to U/Uq, 
is not valid for the small values of Z and used in figs. |] and [7]. This conclusion 
follows from fig. || when the results for U/Uq according to eq. (^) and derived 
from the PA-scheme using eq. ( |Kf ) as static input are compared. Although the 
same approximation is employed for g(r), the two results for U/Uq do not agree, 
since the approximations in truncating the series expansions of the hydrodynamic 
mobility tensors after the terms of (9(r _1 ) and 0(r~ 20 ), respectively, are not the 
same. On the other hand, we have shown in fig. || for highly charged particles 
that the first few terms in the series expansion of eqs. (PJlOl) are sufficient to 
calculate U/Uq. As pointed out before, this is due to the existence of an extended 
correlation hole in the g(r) of highly charged particles in deionized suspensions 
(cf. fig. p]). The correlation hole gives rise to a fast concergence of the integrals in 
eq. ([LTD. However, in case of the weakly charged particles considered here, there 
is no correlation hole present, as can be seen, e.g., from fig. [7]. Therefore, a large 
number of terms in the series expansions of eqs. (PJl^) are needed for weakly 
charged particles to calculate U/Uq. In fact, as discussed before in the case of 
hard spheres, the terms of 0(r~ 3 ) and 0(r~ 4 ) in the series expansions in eqs. 
(|9|,|10D contribute in approximately the same weight to U/Uq, still remarkably 
large compared the the leading Oseen-term of (9(r _1 ), considered only in the 
electrostatic contribution in the PD-result eq. (|20|). 

To summarize the last presented results, we have shown that the range of 
validity of eq. ( PU| ) is restricted to values of <fi considerably smaller than 10~ 3 , 
where the macroion radial distribution function is described very accurately by its 
linearized zero-density form. Furthermore, for the PD-result eq. (^) to be hold, 
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the effective charge Z has to be small enough that the approximation in eq. (O) 
can be used to calculate U/U , but still not so small that the actual g{r) is too 
close to its zero-charge limiting form go{r) = 0(r — 2a). Otherwise, the omission 
of higher order hydrodynamic terms of 0(r~ 3 ) to the electrostatic contribution 
in eq. (^) ceases to be a good approximation. These severe restrictions limit the 
range of validity of eq. (^) to values of the system parameters most probably 
not accessible in sedimentation experiments. Moreover, the relative differences 
in U and Uq become small in the parameter range where eq. (p0|) should apply. 

Let us now turn to some remarks on the 1/,2 -scaling of U/Uq — 1 proposed 
in eq. (^l|). Since eq. (|2T|) is derived from the PD-result eq. (p0[), the be- 
fore discussed restricted range of validity of eq. (|20|) applies also to eq. fl2ip . 
Furthermore, in eq. (pi]), the terms linear in are neglected against the term 
proportional to 1//2 . We have shown in fig. |6]that this additional approximation 
becomes invalid already at small volume fractions, typically ~ 10~ 3 . For larger 
0, the two results for U/Uq in eq. fl20| ) and eq. (^) start to deviate strongly, since 
the contributions to U/Uq linear in (especially the hard core contribution 6.550) 
are no longer small compared to the term proportional to 1//2 . We therefore also 
expect the 1//2 -scaling of U/Uq — 1 according to eq. (^T|) not to be measurable 
in a sedimentation experiment. Nevertheless, we wish to point out that there is 
a specific range of particle charges Z where U/Uq — 1 indeed scales as 1 / 2 for 
a broad range of volume fractions. This range of Z-values is determined when 
the effective particle charge is subsequently lowered from large values Z > 100, 
where U/Uq — 1 scales like 1 / 3 , to Z = 0, where U/Uq — 1 behaves linearly in 
according to eq. (p]). To demonstate the occurance of a 1//2 -dependence of 
U/Uq — 1, we have redrawn in fig. [S] the RMS A- PA result from fig. [5] for Z = 5 
together with the corresponding results for Z = 6 and Z = 7. Obviously, all 
three curves are well described by the form U/Uq = 1 — p<f) a with a parameter a 
very close to 1/2. Such a 1 / 2 -dependence of U/Uq might indeed be measurable in 
sedimentation experiments on dispersions with suitably chosen particle charges. 
However, we stress that the occurance of the parameter a close to 1/2 for a broad 
range of volume fractions as dipicted in fig. |8] is completely different from the 
prediction eq. which has been shown to be valid only for very small 0. 

Finally, we shortly comment on the dependence of our results on the system 
parameters T and e, which are held fixed in our calculations. Obviously, these 
two parameters enter into our calculations mainly by determining the Bjerrum 
length Lb, which has a rather large value Lb = 26.12nm due to the small dielec- 
tric constant e = 2.183 used in this work. As consequence, also the ratio Lb /(2a) 
is comparatively large. This ratio determines together with Z the strength of the 
electrostatic repulsion between the particles. For that reason, one might argue 
that our dicussion concerning the use of eq. (|T9|) for g{r) in calculating U/Uq is 
not longer valid for smaller values of Lb- However, explicit PA-scheme calcula- 
tions for system parameters as used in Ref. [16[ and/or a simple estimate show 
that our conclusions derived above remain valid even for considerably smaller 
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values of Lb'- Consider a value of Lb ten times smaller than the one used in 
the above presented calculations, e.g. a value like in sedimentation experiments 
on charged silica particles in ethanol fTEfl . Then, according to eqs. (|3|,|4]), one 
might choose a particle charge approximately y/lO ~ 3.16 times larger than in 
our calculations to achieve the same strength of repulsion between the particles. 
As shown above for our choice of system parameters, the use of eq. (|T9|) for g(r) 
is a poor approximation even for low volume fractions <fi < 0.01 and a particle 
charge Z = 2. Therefore, even for a system with Lb ten times larger, eqs. ( pO| , pT| ) 
are applicable only for particle charges considerably smaller than Z = 6. This 
slightly extended parameter range for eqs. (pO| , pT|) to be valid is most probably 
still too restricted to be experimentally accessible. 

5 Concluding remarks 

We have presented theoretical results for the reduced sedimentation velocity U/Uq 
of monodisperse charged suspensions in dependence of the volume fraction and 
of the particle charge number Z. Our theoretical model for U/Uq is based on the 
effective macroion fluid model and on the assumption that pairwise additive HI 
prevails at sufficiently small <fi. The numerical results for U/Uq in case of strongly 
charged particles at low salinity are well parametrized by the form 1 — p0 1//3 . 

We have shown that the exponent 1/3 and the value of the charge- independent 
parameter p — 1.8 can be quantitatively understood in terms of a model of effec- 
tive hard spheres of radius aEHS which depends on the volume fraction. Using 
Percus-Yevick input for the static pair correlation function of hard spheres, the 
EHS model can be treated analytically and leads to a value of p very close to 1.8. 

It was further demonstrated that U/Uq — 1 can scale like 1//2 in case of 
dilute suspensions of very weakly charged particles. This peculiar volume fraction 
dependence of U/Uq derives from an expression given by Petsev and Denkov 
[[IT],^] , when it is further assumed that the ionic strength in the system is mainly 
due to counterions. Our numerical calculations reveal that both the original 
expression for U/Uq given by Petsev and Denkov, and the derived expression 
showing the 1 / 2 -scaling of U/Uq — 1, are only valid for very small particle volume 
fractions, which we expect not to be accessible in conventional sedimentation 
experiments. 
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Figure captions 



Figure 1: RMS A radial distribution function g(r) (dashed line) for a suspension 
of charged particles with system parameters a = 695 nm, e = 2.183, T = 293 
K, n s = (i.e. no added salt), Z = 150, and (f> = 0.08. Also shown is the 
radial distribution function of the EHS model in PY approximation (solid line). 
The effective volume fraction is <Pehs — tt/6, corresponding to the identification 
^■Ojehs — r m = r (see text). The inset shows the peak position r m of the RMSA- 
g(r) for the same system paramters, normalized by the mean particle distance r, 
as function of 0. 



Figure 2: PA-scheme result for the reduced sedimentation velocity U/Uq of 
strongly charged particles in a deionized solvent (solid line). The result is well 
parametrized by the form 1 — 1.80 34 . System parameters as in fig. |l|. For 
comparison, the dashed-dotted line represents the lowest-order density result for 
hard spheres according to eq. ([!]). 



Figure 3: PA-scheme results for U/Uq for various values of the effective charge 
number Z as indicated in the figure. All other system parameters as in fig. |l[ 
Notice that U/Uq is nearly independent of Z when Z > 100. 



Figure 4: Dependence of U/Uq on various two-body contributions to the hydro- 
dynamic mobilities. Included terms of the series expansions in eqs. (PJTOD and 
used in the PA-scheme as indicated in the figure. System parameters as in fig. ^| 
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Figure 5: U/Uq versus <ft for particle charge Z = 5. Other parameters chosen as in 
fig. Shown are the PA-scheme results obtained using as static input the RMSA- 
g(r) (solid line), the zero-density form g(r) = 0(r — 2a) exp (—/3u e i(r)) (dashed 
line), and the linearized zero-density form of g(r) according to (|P9|) (dashed- 
dotted line). Further displayed is U/Uq corresponding to eq. ( P0"D (dotted line). 



Figure 6: U/Uq versus <fi for particle charge Z = 2. All other parameters as 
in fig. [I]. Shown are PA-scheme results obtained using as input the RMSA-g(r) 
(solid line), and the zero-density form g(r) = 0(r — 2a) exp (— (3u e i(r)) (dashed 
line). The latter is practically identical with the PA-U/Uq calculated with the 
linearized g(r) of eq. (0). Further displayed are results for U/Uq according to 
eq. ( PU[ ) (dashed-dotted line) and eq. ( pij ) (dotted line). The thin solid line 
represents the hard sphere result of eq. (p. 



Figure 7: Radial distribution functions g(r) for Z = 2 and <fi = 0.005, obtained 
in the RMSA (solid line) and given by g{r) = exp(— (3u(r)) (dashed line). Pa- 
rameters aside from Z and as in fig. [I]. The inset shows the corresponding 
functions r(l — g(r)). 



Figure 8: RMSA-PA-scheme results for U/Uq at particle charges Z = 5, Z = 6, 
and Z = 7. All other parameters as in fig. [I]. The result for Z = 5 (solid line) 
is well parametrized by the form U/Uq = 1 — 2.180 ' 56 , whereas the results for 
Z = 6 (dashed line) and Z = 7 (dashed-dotted line) are very well fitted within the 
shown volume fraction range by U/Uq = 1 — 2.O60 0,53 and U/Uq = 1 — 1.970 0,50 , 
respectively. 
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